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We consider averaged shelling and coordination numbers of aperiodic tilings. 
Shelling numbers count the vertices on radial shells around a vertex. Coordina- 
tion numbers, in turn, count the vertices on coordination shells of a vertex, defined 
, via the graph distance given by the tiling. For the Ammann-Beenker tiling, we 

On ■ find that coordination shells consist of complete shelling orbits, which enables us to 

, calculate averaged coordination numbers for rather large distances explicitly. The 

■ relation to topological invariants of tilings is briefly discussed. 

:s 

Key words: Aperiodic order; Shelling numbers, Coordination numbers; Averages 
1. Introduction 

^ . Many combinatorial questions from lattice theory are best extended to aperiodic system 

^ I by using an additional averaging process. In particular, this is the case for the shelling 

problem, where one asks for the number of vertices on spherical (circular) shells. In 
this contribution, we consider an extension of the shelling problem to the setting of 
more general distances, and give examples for the coordination number case [HI El) which 
corresponds to the graph distance in a tiling. 

To keep things simple, we explain our approach for cyclotomic model sets in the 
Euclidean plane ~ C, with co-dimension 2 (referring to the so-called internal space). 
Here, following the algebraic setting of Pleasants ^Hj, one starts from a set of cyclotomic 
integers, L = Z[^„] with ^„ = e^'^*/" and suitable n, which is the set of all integer linear 
combinations of the regular n-star of unit length. Note that one can choose a Z-basis of 
0(n) elements, where (j) is Euler's totient function [2]. 

This setting is equipped with a natural *-map, defined by a suitable algebraic conju- 
gation (such as ^5 1-^ 1-^ ^g, and ^ C12 examples discussed below, together 
with the canonical extension to all elements of L). The set L := {(x, x*) \ x £ L} is then 
a lattice in M'^^"), the so-called Minkowski embedding of L. A model set A is now a 
set of the form 

A = {x e L\ x*" e Q}, (1) 
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Figure 1: Patches of the Ammann-Beenker tihng (left) and the shield tiling (right). 



or any translate of it, where the window is a relatively compact subset of internal 
space with non-empty interior. A natural choice for 17 that preserves n-fold symmetry is 
a regular n-gon, which leads to regular model sets (i.e., the boundary df2 has Lebesgue 
measure 0). More precisely, we focus on the generic case {df2 f] L* = 0), where A is 
repetitive and its Ll-class (consisting of all locally indistinguishable patterns) defines a 
uniquely ergodic dynamical system [ITj. This includes the examples of figure ^ 

Let d{x, y) denote any translation invariant distance between x and y. Due to unique 
ergodicity, combined with finite local complexity, we know jl7l |Sj that the averaged 
number s^{r) of vertices on a d-shell of radius r is determined by a sum over a finite 
number of patches, weighted by their frequencies, which exist uniformly. Since we work 
with model sets, this can be further reduced to a sum that involves only admissible pairs 
of vertices with the correct distance. The frequency of a pair (x, y) with x,y G L is given 
by the scale-independent autocorrelation coefficient ^{x — y), where [S] 



with 1^ denoting the characteristic function of the window. The normalisation is such 
that z^(0) = 1, i.e., we count per point of A rather than per unit volume. Clearly, 
i/{z) = v{—z), and further identities may occur as a result of the symmetries of the 
window. In general, the averaged number reads 



which can then be further simplified by means of a standard orbit analysis 

In view of this derivation, it is reasonable to consider a d-shell to be the collection of 
all pairs (x,y) of a given distance d{x,y) together with their frequencies z^(x — y). 

The coordination problem is now the extension of the shelling problem to a different 
distance concept, based on the graph distance in the tiling under consideration. The 
graph distance of two vertices is defined as the minimum number of edges in a path 
linking the two vertices. For simplicity, we restrict our discussion to examples with a 
single edge type (i.e., all edges have the same length), though various extensions are 
possible. 
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Table 1: Averaged coordination numbers of the Ammann-Beenker tiling. 
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There is an important connection between averaged coordination and shelling num- 
bers, which stems from the relation between graph and Euclidean distances. On a given 
coordination shell, the vertices appear in symmetry orbits of the underlying tiling; these 
vertices are orbitwise distributed over finitely many Euclidean shells. Within one orbit, 
every vertex contributes the same amount to the averaged numbers. Conversely, the 
vertices on a Euclidean shell orbitwise belong to (possibly different) coordination shells. 
Here, the term orbit simply refers to the orbit of a point under the point symmetry group 
of the pattern under consideration. 

2. Examples 

In what follows, we concentrate on three particular examples, with 8-, 10-, and 12-fold 
symmetry. For n = 8, we consider the Ammann-Beenker tiling, obtained by the above 
construction with a regular octagon of edge length 1 as the window. The special role 
of \/2 reflects the fact that Z[-v/2] = Mn ^[^g], see [S] for details. This reference also 
contains the known results on the circular shelling. 

For this tiling, the connection between shelling and coordination numbers is rather 
advantageous, because coordination shells comprise only complete circular shells. This 
is a consequence of the fact that the four directions of the edges in the tiling form a 
Z-basis of the underlying module ^[^g]) which is possible because <^(n) = n/2 for n = 8. 
This means that for a given distance vector, the number of steps along each direction 
is uniquely determined. Moreover, there always exists at least one path along edges of 
the actual tiling that is admissible (in the sense that never has to 'backtrack' along the 
path). This is a higher-dimensional analogue of the corresponding (trivial) situation for 




Figure 2: The first 400 coordination numbers for the Ammann-Beenker tiling (left) and 
the fluctuations in the differences of successive coordination numbers (right). 



the silver mean chain (when viewed as a cut and project set obtained from a rectangular 
lattice), which fails for other tilings. 

The averaged coordination numbers can therefore be calculated by identifying the 
contributing circular shells, and summing the corresponding averaged shelling numbers, 
which can be obtained as described in 0. Following this approach, we calculated the 
first few hundred averaged coordination numbers for the Ammann-Beenker tiling, see 
tableland figure 121 

As for periodic planar lattices, the averaged coordination numbers Sc{k) grow, on 
the average, linearly with the number of steps k. However, a closer inspection of the left 
part of figure 121 shows that the growth rate fluctuates, and that the data points do not 
lie on a single line, but inside a sector bounded by two lines of slightly different slopes. 
An even closer inspection reveals that the fluctuations of Sc{k) follow a sophisticated 
pattern, as displayed in the right part of figure (21 which shows the differences Asc{k) = 
Sc(A;-|- 1) — Sc{k) of consecutive averaged shelling numbers. The resulting pattern appears 
to comprise a number of sinusoidally varying curves. This might be caused by the 
variations of overlap areas of the window with shifted copies of itself, which enters the 
computation of patch frequencies; a closer investigation of this phenomenon might lead 
to interesting results. 

A construction of the rhombic Penrose tiling as a cyclotomic model set with four 
components is described in [J]. Due to the necessity of four windows (which are also 
present in the non-minimal embedding via the lattice Z^), the determination of averaged 
quantities is technically more involved. Table [21 recalls some results obtained earlier in 
1^1. In this case, ^[^5] nM = Z[r], where t = {1 + \/5)/2 is the golden ratio, which is the 
relevant algebraic integer here. 



Table 2: Averaged coordination numbers of the rhombic Penrose tiling. 
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Table 3: Averaged coordination numbers of the twelvefold symmetric shield tiling, split 
into the contributions from circular shells of radius r. 
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As an example with 12-fold symmetry, we consider the so-called shield tiling 
of figure n It is obtained from Z[^i2] choosing a regular dodecagon of edge length 1 
as window, see 4 for details. In this case, Z[Ci2] n M = Z[^/3]■, examples are listed in 
table El Note that a single circular shell can contribute to several coordination shells. 

3. Frequency modules 

It is remarkable that the averaged coordination numbers are special algebraic integers in 
all three examples. From the cut and project method, in conjunction with equations ((21) 
and ©, it is clear that these numbers must be rational, i.e., elements of the corresponding 
cyclotomic field Q(Cri). This follows from the computability of the integrals in equation 
((JJ within these number fields. 

The further restriction to algebraic integers is due to a special structure of the fre- 
quency module of the tiling, i.e., the Z-span of the frequencies of all finite patches in 
the tiling. Since our averaged quantities are simple integer linear combinations of patch 
frequencies, a 'quantisation' of the latter to integers implies the result for the former. 
This phenomenon has been observed before several times, and for various related prob- 
lems [SI El 121 El. The proof relies on the topological structure of the compact Ll-class, 
viewed as a dynamical system under the translation and/or inflation action [H l8l (Toj. 
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